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ABSTRACT 

A thin dark thread is observed in a UV/EUV solar jet in the 171 A, 193 A and 211 A and 
partially in 304 A. The dark thread appears to originate in the chromosphere but its temperature 
does not appear to lie within the passbands of the Atmospheric Imaging Assembly (AIA) onboard 
the Solar Dynamics Observatory (SDO). We therefore implement solar magnetoseismology to 
estimate the plasma parameters of the dark thread. A propagating fast kink (transverse) wave 
is observed to travel along the dark thread. The wave is tracked over a range of ~ 7000 km 
by placing multiple slits along the axis of the dark thread. The phase speed and amplitude of 
the wave are estimated and magnetoseismological theory is employed to determine the plasma 
parameters. We are able to estimate the plasma temperature, density gradient, magnetic field 
gradient and sub-resolution expansion of the dark thread. The dark thread is found to be cool, 
T < 3 X 10^, with both strong density and magnetic field gradients. The expansion of the fiux 
tube along its length is ~ 300 — 400 km. 

Subject headings: magnetohydrodynamics (MHD) — plasmas — magnetic fields — Sun: oscillations 



1. Introduction 

Short-lived, transient phenomena are common- 
place in the solar atmosphere. In particular, large 
jet like features are commonly seen to erupt from 
the lower to mid solar atmosphere, ejecting ma- 
terial into the corona. These phenomena are 
grouped under the term solar jets, which covers a 
wide range of explosive events that have been ob- 
served over the last century, includ ing Ha surges 



ported a number of times among the jets seen in 
UV/EUV is a large-scale torsional/helical mo- 
tion fe.g. fShimoio et al I Il996l: IPike Mason 1 



(iNewton I Il934f ). macro-spicules (iBohlin et al 



1975h. UV/EUV (iBrueckner fc Bartoe I ll983h. X 



(jPike Mason Ill998h . 



rav jets ([Shibata et al. 1992) and solar tornadoes 



19981: IPa tsourak os et al. I l2008l : IlTu et al. I hom 
Kamio e t al. 20l3). Proposed reconnection mod- 
els suggest that closed twisted field lines re- 
connect with surrounding open field lines (e.g. 
Shibata fc Uchida Ill985l : IPariat et al. Il2009[ ). 

Whatever the exact excitation mechanism of 
these jets is, it appears to coincide with the ex- 
citation of a number of different types of wave 
phenomena. Firstly, there is the large-scale tor- 
sional /helic al motion reported in both observa- 



It is believed that these events are driven 
by relatively localized reconnection events (e.jg; . 
Yokovama fc Shibata llT995l . Il996l : lHarrisonllT999[ ). 
where emerging magnetic fiux reconnects with 
existing fields. One feature that has been re- 



tions fe.g. iKamio et ah] 2010[ ) and simulations 
(e.g. IPariat et al. I l20Q9f ). whose interpretation 
is still in debate (t he candidates being torsional 



Alfven waves (e.g. IPariat et alH l2009[ ). or heli- 
cal kink waves, M.S. Ruder man - private com- 
munication). The second type of reported wave 



1 



phenomena is the t r ansverse (fast ki nk) wave 



([Cirtain et aTl I2QQ7I : iLiu et al. I |2QQ9[ ). How- 



ever, in the majority of observations of trans- 
verse (fast kink) waves in jets, the wave phenom- 
ena, in our view, was wrongly interpreted as an 
Alfven wave, whose features are quite distin ct 
from the transverse kink wave (lErdelvi FeduiTI 



2007 : Van Doorsselaere et al. ]l2008'). The trans- 
verse kink waves are normally reported in a bright, 
fast moving part of the jet which is common to 
the X-ray jets. 

The study of wave phenomena in the solar at- 
mosphere proves useful as demonstrated by the 
new and expanding field of solar magnetoseis- 
mology, which exploits observed wave phenomena 
to determine hard to measure or unmeasurabl e 
plasma parameters (Uch ida 1970l : Roberts et al 



Ll984; Erdelvi 2006a. b). Recent advances in ob- 
serving technology has demonstrated the ubiqui- 



e.g. Baneriee et al. 


20071: Andries et al. 


2009; 


De Moortel 1 20091: 


Ruderman & Erdelvi 


2009; 


Tarovan & Erdelvi 


20091 for reviews of waves 



in the solar atmosphere). However, magneto- 
seismology has largely been confined to observa- 
tions of structures with relativel y long li fe times 
e.g. coronal loops, (e.g. Tarovan fc Bra dshaw 



2QQ8l:IVerth et al. kOOSHTaroyan fc Erdelvi 11200^ 
Morton fc Erdel vi I [2009ai I2OIOI : IVerth et alT 



2010; Terradaseta 



i 



201 ll). promin ences (e.g. 



Pinter et al. I I2OO8I : ISoler et al. 1 12010[ ). although 



recent papers have extended solar magnetoseis- 
mology to transient events such as X-ray jets 
(^Vasheghani Farahani et al. 1120091), coronal waves 



West et al. I2OIII) and spicules ( Zaqarashvili fc Erdelyi 
[2009; Verth etalJ^Oll). 

Why should magnetoseismology be used to 
study solar jets? One reason is that jets are 
short lived events with even shorter lived com- 
ponents, e.g. a ssociated spicu l e-like features at 
the footpoints (Sterling et al. I l2010h . Regular 
methods of determining plasma properties, e.g., 
magnetic field gradients by spectropolarimetry, 
can requi re fairly long in tegration times. For 
example, ICenteno et al. I ( [201 0) attempted to 
measure spicule magnetic field strengths. How- 
ever, large signal to noise ratios in off-limb 
measurements meant integration times of 45 
minutes. Spicules exist only for 5-15 minutes 
(| Zaqarashvili fc Erdelyi 1 120091 ) so such measure- 



ments are much longer than the lifetimes of the 
objects they are trying to study. Another rea- 
son is that properties of dark features, i.e., little 
or no intensity at certain wavelengths, cannot be 
measured from the spectral information at those 
wavelengths. Solar magnetoseismology negates 
these difficulties as the changing properties of the 
wave, i.e. period, phase speed, amplitude can 
reveal the required information on plasma strat- 
ification properties, i.e. magnetic field gradient, 
density scale heights, etc. Thi s approach has been 
demonstrated successfully bv lVerth et J~l (l201lh 
in the study of spi cules. Further to th is, solar jets 
are multi-thermal ( Kamio et al. iBoiOf ). X-ray jets 
in particular have a rapid, hot ejection of mate- 
rial seen in X-rays and a much cooler ejection of 
material that m akes up the jet obse rved in EUV 



lines (see, e.g.. iKamio et aLl l2010[ ). Obtaining 



values of plasma parameters, e.g., temperature, 
and variation of plasma parameters with height, 
e.g., density, magnetic field, of multi-temperature 
features from the spectroscopic data available in 
the narrowband filters, such as those onboard the 
Solar Dynamics Observatory (SDO), proves diffi- 
cult if not impossible. Solar magnetoseismology, 
again, negates this problem. For example, magne- 
toseismology can be used to determine the density 
and magnetic field gradients, and temperature, in 
theory, of a feature that appears dark in an SDO 
filter. Little or no information could be obtained 
about the plasma in the object of interest from 
the filter response itself. Finally, magnetoseismol- 
ogy allows the determination of fine structure at 
sub-resolution scales. It is at these small scales 
that wave heating may occur, so the study of such 
length scales is vitally important. 

In the following we present the observations of 
a propagating kink wave observed in a structure 
associated with a solar (UV/EUV) jet. The struc- 
ture is a dark inclusion, in a relatively hot jet, 
whose presence appears not to have been seen be- 
fore (at least in EUV lines). We then use magneto- 
seismology to estimate the temperature, the den- 
sity stratification and expansion of the observed 
magnetic structure. We compare the obtained 
results with the measured features of the larger 
jet structure, magnetic field extrapolations for the 
quite Sun magnetic field and chromospheric esti- 
mates. 
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Fig. 1. — First and second rows: Images showing the morphology of a large jet (observed on 20 January 
2011) over its life time taken by SDO/AIA. The panels in the top row are from the 304 A channel taken at 
09:16, 09:18, 09:25 and 09:33. The images on the bottom row are from the 171 A channel taken at the same 
times. The bright loop is seen in both channels at 09:16, suggestive of reconnection. A dark feature is seen 
at 09:18 at the same position as the bright loop appeared. This dark feature grows into the dark thread seen 
in 171 A at 09:25. Third and fourth rows: Context images of the jet from STEREO. The first three panels 
in each row are from STEREO A and the last is from STEREO B. The top row are images from the EUVI 
304 A channel and are taken at 09:06, 09:26, 09:36 and 09:36. The bottom row are images from the EUVI 
171 A channel and are taken at 09:03, 09:28, 09:33 and 09:33. 
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2. Observations 



2.1. Solar jet 



The observations began at 09:00 UT on the 20 
January 2011 and last an hour and a quarter tih 
10:15 UT, on the south west hmb us ing the Atmo- 
spher ic Imaging Assembly (AIA) (|Lemen et al 



20111 ) onboard SDO which has a spatial resolu- 
tion of ~ 0.6 arcsec per pixel or two-pixel spatial 
resolution of ~ 870 km and a cadence of ~ 12 s. 
The solar jet was observed mainly in the 304 A 
spectral line with faint emission signatures in 171 
A, 193 A and 211 A, close to the southern polar 
coronal hole. The time series were obtained from 
the SSW cutout service which provides level 1.5 
data that has already been corrected for flat-field, 
despiked and the images co-aligned. The line in 
which the jet appears brightest is the 304 A which 
has a peak temperature of ~ 10^*^ K (see first 
row in Figure [T]). Our main interest lies in a dark 
thread in the 171 A line (second row in Figure [T]); 
however, we will also describe the part of the jet 
seen in the 304 A channel as it will be important 
later. 

In Figure[T]we also provide images from the Sun 
Earth Connection Coronal and Heliospheric In- 
vestigation (SECCHI) telescopes onboard the So- 
lar Terrestrial Relatio ns Observatory (STEREO) 
( Howard et al. |[2008[ ) to show the orientation of 
the jet in the line of sight with respect to SDO 
(third and fourth rows). The bright UV/EUV 
portion of the jet has a 18 degree angle with 
the normal to the Sun's surface; however, it is 
not possible to determine the angle between the 
dark thread and this normal. If it is similar to the 
UV/EUV jet then any projection effects should 
may be relatively small but should be kept in 
mind. There is also a clear emission in the 171 A 
filter from the jet, which appears as a relatively 
distinct spike in the quiet corona. The EUV emis- 
sion is relatively dim suggesting the bulk of the 
material in the jet is less than 1 MK, although 
due to blending of the 171 A with cooler lines the 
actual v alue could be sig!:nificantly less than this 



(see, e.g. lDel Zanna fc Mason Il2003[ ). 



The jet occurs close to the limb at the southern 
pole of the Sun and is situated just outside a coro- 
nal hole in a relatively quiet region. We do not go 
into a detailed discussion here of how the jet may 
be driven; however, one possible scenario is the jet 



could be generated by a reconnection event in the 
lower solar atmosphere. This will be the subject 
of a further investigation. In Figure [T] we show the 
morphology of the jet over its lifetime at various 
times. 

In the 304 A line we see a small bright colli- 
mated jet rise from the surface. The jet expands 
with height and time and during its main phase, 
highlights a funnel shape. We place cross-cuts of 
70 pixels in length and one pixel in width per- 
pendicular to the central axis of the jet at heights 
of ^ 10 Mm, ^ 14.3 Mm and ^ 15.7 Mm above 
the limb and use linear interpolation to obtain an 
X — t data set. We then measure the width ex- 
pansion of the jet over time at various heights as 
a function of time. Figure [2^. (N.B. in all figures 
from now on the times are given in seconds from 
9:22 UT, which is when the kink wave is first ob- 
served. The benefit of this is clear in Section 4 on 
magnetoseismology.) The jet is clearly defined due 
to the jet material having a much greater intensity 
than the surrounding chromospheric material and 
the jet edges are determined by eye. We suggest 
there could be an error of < 10% on these mea- 
surements. Figure [2^ shows that there is a sim- 
ilar trend for the expansion of the jet over time 
at all heights. As certain reconnection scenarios 
involving reconnection between emerging and ex- 
isting field describe an untwisting magnetic field, 
it is possible that the expansion of the jet cross- 
section in time is partly due to to the untwisting 
of the field. The expansion of the fiux tube in Fig- 
ure [2^ is not constant but goes through a number 
of phases. 

At around 16 Mm above the surface the jet 
experiences a whip-like effect (see F i gure [J) sim- 
ilar to that reported in iLiu et al. I (|2009[ ). The 



plasma above the whip appears colHmated. Once 
the plasma has finished being ejected into the at- 
mosphere, a portion of it returns to the surface. 
This falling material highlights the magnetic field 
and it follows a collimated path back to the chro- 
mosphere suggesting the complete untwisting of 
the magnetic field has occurred. The funnel shape 
of the magnetic field has now evolved into a col- 
umn. 

Further to this, we calculate the median value of 
intensity fiux, F, in the 304 A channel of the pix- 
els in the cross-cut placed at 10 Mm, which were 
determined to lie within the jet boundaries. This 
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is plotted as a function of time in Figure [2]3. We 
do this because the intensity is dependent upon 
thermodynamic quantities, e.g., temperature and 
density, so provides an estimate of how the den- 
sity/temperature in the jet evolves as a function 
of time. 

In the 171 A channel the jet appears completely 
different to that seen in 304 A. At 09:16 a bright 
loop is clearly seen with the eruption starting two 
minutes after. A thin, dark thread, with a width 
at resolution scales, is then seen to be ejected into 
the atmosphere along with a very small amount 
of bright material. The dark material appears to 
come from the chromosphere, appearing (appar- 
ently) from behind the bright loop (09:18) (Fig- 
ure d]) . This is the first appearance of the dark 
thread that is shown in the image at 09:25. The 
dark thread is co-spatial and co-temporal with 
the UV/EUV component of the jet in the 304 A 
channel but is barely visible amongst the hotter 
plasma. In the middle right panel of the 304 A 
images a bright, thin thread of plasma is seen to 
highlight the right hand jet edge. This can just be 
made out in the corresponding 171 A image run- 
ning along the righthand side of the dark thread. 
Towards the base of the bright thread in the 304 A, 
a small dark inclusion can be seen which is the 
dark thread (e.g. third panel, top row. Figure [T]). 
This suggests that the thin, dark thread is a flux 
tube (or flux tubes) contained within a large col- 
lection of flux tubes, i.e. the jet seen in the 304 A 
channel. The dark thread survives for around 10 
minutes before the material either traveling higher 
into the atmosphere or returning to the surface. 
Further, the dark thread jet has an angle of ~ 35 
degrees to the normal from the solar surface. The 
plane in which this angle is measured, i.e., the 
plane of view of SDO/AIA, is almost perpendicu- 
lar to the plane of view of STEREO. Hence, the 
angle measured for the bright part of the jet in 
STEREO is measured in a different plane to this 
angle. 

Aside from the dark thread, only a small 
amount of bright plasma is seen to be ejected 
into the atmosphere in 171 A. This suggests that 
the density and/or temperature of the dark thread 
plasma must be different from the material seen in 
the 304 A channel. Due to its appearance as a dark 
inclusion we suggest the material is much cooler 
and/or denser (hence, optically thicker) than the 




other jet material. The peak sensitivity of the 
304 A channel is T ~ 10"^-^ K, so the material in 
the dark thread must be at least T < 5 x 10^ K. 
This, along with the appearance of the bright loop 
at the jet footpoint, supports the popular conjec- 
ture that the reconnection event responsible for 
the jet occurred at chromospheric heights. The 
dark thread could be si milar to the cool chromo- 
spheric jet observed by iLiu et aL ] (2009). One 
dif ference between t he jet observed here and that 
(120091) is that the ejected material 
(l2009h returns to the surface along a 
similar path to that of its ejection. Here, some of 
the material is ejected higher into the atmosphere 
while the rest returns to the surface along a dif- 
ferent path at least a few Mm from the ejection 
site. 

The dark thread has a width close to the reso- 
lution limit so the expansion that is measured for 
the UV/EUV jet structure is not seen in the dark 
thread. The outflow of the dark thread material 
is also a lot slower than outflow of bright material 
seen in the 171 A. In 171 A we can measure the 
speed of the outflow for the dark thread by using 
two methods. For the initial rise of the dark fea- 
ture we use a running difference movie to track the 
dark thread front and obtain average rise speeds 
of 63 kms~^. After the initial rise we employ a 
method of contouring features of different inten- 
sity and tracking th e intensity features thro ugh 



2001L for 



the images (see, e.g. IWinebarger et al. 
details). The flow in the dark thread is given as 
an average speed along the thread of the different 
intensity features identified and tracked through 
the thread. Here, we term the initial rise or possi- 
ble 'extension'of the dark thread as outflow. In the 
following section we need to subtract any apparent 
motion, both flow and extension, from measured 
phase speeds. We use the term outflow to describe 
the observed upwards motion of plasma, making 
no definite distinction between possible extension 
or flow. 

In Figure [3] we show the outflow speeds along 
the dark thread as a function of time. The speed of 
the outflow hes within the ran ge of outflow speeds 



reported in iLiu et al. I (|2009l ). We then smooth 
the measured values of flow with a 5-point boxcar 
average and fit a quadratic profile to the smoothed 
points which is the dash dot Hue in Figure [3l It 
would appear that there are two distinct phases to 
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Fig. 2. — (a) Expansion of the jet width in 304 A as a function of time plotted for different heights along the 
jet. The solid line is taken at 11 Mm from the base, the dashed line 14 Mm and the dotted line 15 Mm. The 
time is given in seconds from 09 : 22 UT. (b) Average intensity of the jet in a cross-cut taken at a height of 
10 Mm above the solar surface as a function of time in 304 A. The flux is measured in terms of the number 
of photon counts per second. 



the flow; a phase with decreasing flow speed and 
a phase of almost constant flow speed. Separat- 
ing these into two different phases we construct 
quadratic fits to the two phases (solid lines). The 
second phase can be seen to have an almost con- 
stant flow speed. It is unknown why the outflow 
rate appears to have these two distinct phases, but 
is probably connected to the nature of the explo- 
sive (reconnection) event. 

2.2. Transverse wave 

We can determine very little about the dark 
thread from the SDO filters due to its low in- 
tensity in all SDO filters. However, we observe 
a propagating transverse wave traveling along the 
thread. The natural interpretation for the wave is 
a propagating transverse wave (propagating kink 
wave). A number of previous authors have re- 
ported observed displacements of the jet axis as 
the Alfven wave. However, Alfven waves are oscil- 
lations of constant magnetic surfaces and do not 
displace the flux tube axis, the only mode that 
causes a visible displacement of the flux tube axis 
is the fast kink mode (see, e.g; . lErdelvi fc Fedun 
2007 : Van Doorsselaere et al. I l2008l for discus- 



sion). The transverse wave should allow us to 
estimate certain plasma parameters by means of 
solar magnetoseismology. First, in this subsection 
we describe how we obtain the information about 
the wave and in the following sections we discuss 
the magnetoseismology. 

We determine the line parallel to the final axis 
of the dark thread and then place 9 cross-cuts 
(linearly interpolated) perpendicular to the line, 
with the cross-cuts separated by 0.97 Mm (see. 
Figure S]). The wave appears at approximately 
half way along the thread, which is approximately 
^ 7 Mm above the limb in 171 A. We do not 
consider projection effects along the line of sight 
for estimates on distances and heights. The angle 
between the normal to the solar surface and the 
UV/EUV jet is approximately ~ 18 degrees (see. 
Figure [T]) , so the correction will only be relatively 
small, ~ 5 % on measured distances. The prop- 
agating kink wave can be seen up to a height of 
~ 15 Mm in the atmosphere. This corresponds to 
the wave traveling a distance of ~ 11.6 Mm along 
the dark thread. An example of the oscillation ob- 
served in one of the cross-cuts is shown in an x-t 
plot in Figure [5^. 
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Fig. 5. — (a) A sample cross-cut x-t plot (cross-cut No. 3). The transversal kink motion of the dark thread is 
clearly seen, (b) and (c) An example of the fitting procedure for the third cross-cut. (b) The middle panel 
shows the data points (crosses) over plotted with a smoothed profile (solid line). The dotted line shows 
the linear fit to the wave, (c) The right hand panel shows the wave with the linear fit subtracted. The 
symbols indicate the first zero-crossing {square)^ the maximum {triangle)^ the second zero crossing {cross) ^ 
the minimum {circle) and the third zero crossing {diamond). 
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Fig. 3. — Measured values of flow as a function 
of time along the dark thread in the 171 A. The 
crosses are the measured values while the solid 
Hues correspond to the quadratic fits to the mea- 
sured values smoothed with 5 point boxcar aver- 
age. The dashed lines are the 95 % confidence 
levels from the fits. The dashed-dot line is a 
quadratic fit to all the data points 
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Fig. 4. — Image of the dark thread showing the 
cross-cut positions. The observed transverse mo- 
tion propagates from cross-cut 1 to 9. 



To obtain information about the wave we deter- 
mine the pixel in each time frame of the x-t plot 
that has the minimum intensity. The dark thread 
has an almost inverse gaussian intensity profile 
across its width and the minimum intensity should 
correspond closely to the minimum of this inverse 
gaussian. We then smooth the data points of the 
oscillation using a 5 point box-car function before 
applying a linear fit to the smoothed data. The 
linear fit is subtracted from the smoothed data to 
remove the trend. An example of these steps is 
demonstrated in Figure [Sb and c. We then deter- 
mine the times at which the zero crossings and the 
maximum and minimum values of amplitudes oc- 
cur in the different cross-cuts. The magnitude of 
the maximum and minimum values of the ampli- 
tudes are also determined. Above the 9th cross- 
cut the strength of the wave signal is weak making 
it difficult to obtain fits to the observed profile. 
We begin by plotting the travel times at which 
the zero crossings and the maximum and mini- 
mum values of amplitudes occur at the different 
heights in Figure [6^. The cross-cut labelled No. 1 
in Figure |4] is the first cross-cut in which we see 
the oscillation; hence, in Figure [6^ we measure 
distance from this cross-cut. The crosses corre- 
spond to the values taken from the observations 
and the Hues correspond to a quadratic fit to the 
data points of each feature. The gradient of the 
quadratic fits gives the speeds at which the zero 
crossings and the maximum and minimum values 
of amplitudes travel along the dark thread and is 
shown in Figure [7^. 

One might expect that the speeds of each point 
should have similar behavior as they propagate to 
the different heights. However, as we have illus- 
trated in Figures [2] and [3l the plasma properties, 
i.e. flow, magnetic field and density, in the larger 
collection of flux tubes (304 A jet) and the dark 
thread are subject to large changes over the life- 
time of the jet/thread. If we assume a linear fit 
so that the phase speeds are constant, we see the 
trend is a decrease in the phase speed over time. 
This appears natural as the phase speed of the 
transverse wave is proportional to the magnetic 
field, so an untwisting field would decrease the 
field strength, hence, the value of phase speed will 
also decrease. N.B. Neglecting inclination of the 
flux tube with respect to the normal to the so- 
lar surface gives errors of ±3 kms~^ which lie well 
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within the errors/confidence levels given by the 
fitting procedure. 

Finally, we plot the size of the maximum and 
minimum amplitude envelopes in Figure [6)3. It 
is seen from the amplitude envelope plot initially 
there is an increase in amplitude with height fol- 
lowed by a decrease. The confidence intervals on 
the fits demonstrate that there is a possibility that 
the amplitude does not necessarily have to be de- 
creasing after 4000 km. It is clear from the con- 
fidence intervals in Figure [GJd that the minimum 
amplitude appears to have large errors associated 
with it. In the following sections we will use only 
the maximum amplitude envelope to increase con- 
fidence in our results. 

3. Theory 

Before we analyze the wave motion any further 
we need to determine which, if any, theory al- 
ready developed is relevant to these observations. 
A complete theory that would describe kink mo- 
tion in the jet scenario has not yet been fully de- 
veloped due to the immense complication asso- 
ciated with including all the topological and dy- 
namic features that the reported jet possess. How- 
ever, we will assume static behavior and adapt the 
theory for a magnetic flux tub e with both mag- 



netic and density stratificatio n (|Verth fc Erdelvi 
boost iRuderman et al. Il2008[ ). If the plasma is 



indeed cool as we suspect, then it may be better 
to u se models inc l uding a partially ionized plasma, 
e.g. JSoler et al. " 



(|2009al lbl) who model prominence 
plasmas. However, the results from these models 
indicate partial ionization, at least in single fluid 
theory, has a small effect on the propagation and 
damping of fast kink oscillations (for typical ob- 
served wavenumbers). 

First, we need to make sure any assumptions 
used in the derivations of the previous theories 
are also adequate for this dark thread: (i) The 
thin tube approximation applies. For this we 
require the wavelength of the oscillation to be 
much greater than the radius. The width of the 
dark thread is measured to be < 500 km. We 
can estimate the wavelength/wavenumber of the 
kink wave. The typical period of the oscillations 
(e.g. see oscillation in Figure [5|) is 360 s. Using 
P = 2t: /uj = \/ Ck and values of phase speed from 
Figure [Zb then, 5400 < A < 18000 km and hence 



R/\ < 0.1. (ii) The plasma- /3 is small. Estimates 
for the sound speed in the chromosphere are in 
the range 5-12 kms-^ (T = 5 x 10^ - 1 x 10^ K) 
while the Alfven speed is between 50 — 200 kms~^. 
The plasma- /3 is given by /3 = 2c^/(7i'^) so esti- 
mates for the jet give j3 < 0.25. From these esti- 
mates we can consider the plasma- /3 in the dark 
thread to be small, (iii) The plasma inside the 
dark thread is much denser than the surrounding 
material, i.e. pi » pe- Although this is not a 
necessary condition, it allows simplification of the 
governing equations and suggests that the proper- 
ties of the internal plasma dominate the wave be- 
havior. The dark thread appears to originate from 
the chromosphere/transition region and ejects the 
material from these regions high into the corona. 
It is generally well accepted that corona is much 
less dense than the lower layers of the solar atmo- 
sphere; hence, pi is at chromospheric/ transit ion 
region densities and pe is at coronal densities. 
However, the plasma at UV/EUV temperatures in 
304 A appears to surround the dark thread, so the 
external plasma could be considered to be at up- 
per chromosphere/transition region temperatures 
and densities. It is then much less clear whether 
pe is less than pi. 

First, the phase speed for a fast kink (or trans- 
verse) mode in a thin tube is given by 



( B,{z,tf + B,{z,tf 



1/2 



(1) 



where the subscripts i and e refer to the internal 
and external values of the plasma quantities. In- 
stead of making any assumptions about the rela- 
tive magnitudes of the densities and magnetic field 
we will rewrite this in terms of the average density 
and magnetic field 



B 



Ck 



(2) 



where p = {pi{z,t) + Pe{z,t))/2 and i? = 
iB,{z,t)^ + B,iz,tf)/2. 
Now, in lRuderman et ah 



(|2008D the kink wave 



equation for a magnetostatic flux tube with both 
magnetic and plasma density stratification in the 
longitudinal direction is given by 



(Pt] 



7? = 0, 



(3) 
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Fig. 6. — (a) Shown are the times at which the zero-crossings and maximum and minimum amphtudes 
appear at the different heights (cross-cuts). The data points (crosses) have been fitted with a quadratic 
profile (solid line) and the 95% confidence interval is shown for the fit. (b) Maximum and minimum values 
of the amplitude envelope plotted as a function of height. The crosses correspond to the measured values, 
the solid line is a quadratic fit to the data points and the dashed lines correspond to a sigma confidence 
interval for fit to the maximum amplitude envelope and 0.5 sigma confidence interval for fit to the minimum 
amplitude envelope. 
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Fig. 7. — (a) Phase speeds for first zero (solid), maximum amplitude (dotted), second zero (diamonds), 
minimum amplitude (dashed), third zero (dashed dot), (b) Measured phase speed minus the measured flow 
for the maximum amplitude (solid line) and second zero crossing (crosses). The dashed lines correspond to 
the sigma value error on the phase speeds. 



where rj = £,±/R{z)^ ^± is the displacement per- 
pendicular to the magnetic field, R{z) is the loop 
radius and uj is the frequency of the transverse 
oscillation. 

First, we introduced a scaled length, ( = ez 
under the assumption the length scale of the 
changes in the background plasma parameters, 
i.e. density, magnetic field, are small compared 
to the wavelength, hence, the kink speed changes 
slowly along a stratified and expanding waveg- 
uide. _\^^jiowa22bL_tiie WKB approximation (see 
e.g., bender fc Orszag in the z direction, 

assuming that 

V = A{C)exp{-ie-^e{C)). (4) 

Substituting Equation @ into (|3]) and taking the 
zeroth order terms, with respect to e, gives the 
frequency 

UJ = Ck{z)k, (5) 

where k = d&/d(. The first order terms provides 
an equation for the amplitude of the wave, namely 



where the dash corresponds to derivative with re- 
spect to C- Solving this and using our definition 



of 1] we obtain 

a = J^m- (7) 

V (jJ 

The frequency is independent of the changes in the 
spatial direction, hence, if we can determine the 
change in phase speed and spatial displacement of 
the thread (amplitude) with height then the ex- 
pansion along the fiux tube can be calculated. 

4. Solar magnetoseismology 

Now, the data we have determined from the 
observations contains a large amount of informa- 
tion. To demonstrate the magnetoseismology here 
we shall only use two of the obtained phase travel 
times, for which the reasons will become clear. 

First, we return to the plot of the phase speeds 
shown in Figure [7^. For a propagating wave in 
a fiowing medium, the measured phase speed will 
be a sum of the kink speed and the fiow speed (if 
the wave is propagating in the same direction as 
the flow). If we subtract the measured value of 
flow from the measured phase speeds then we see 
that the actual phase speeds cover a much smaller 
range of values (e.g. see Figure [Tb). From Equa- 
tion ([T]) we can see the kink speed is also depen- 
dent upon the magnetic fleld and the density. The 
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dynamic changes observed in the large jet could 
also be occurring in the dark thread and could be 
sufficient to explain the remaining differences be- 
tween the phase speeds. 

4.1. Scale height and temperature esti- 
mates using the phase speed 

The UV/EUV jet is clearly highly dynamic, 
where magnetic field and density undergo rapid 
changes. However, if we consider Figure [21 we 
can see that there is a period of time when the 
time- dependent expansion and change in inten- 
sity of the jet is at a minimum {t ~ 170 s to 
t ~ 400 s). In this period of time the total expan- 
sion of the jet is < 5% and the change in intensity 
is < 10%. For this time period, the second zero 
crossing is propagating along the thread. We then 
make the assumption that the activity in the dark 
thread is also at a minimum at this point. If this 
is the case, the change we see in phase speed at 
the different heights is due only to the structuring 
along the thread and not due to dynamic behavior. 
The second zero crossing appears first at around 
156 s (9:24:30 UT), propagates along the cross- 
cuts within 110 s and experiences an increase in 
phase speed as it travels along the thread. Us- 
ing the estimate of flow speed along the jet during 
this period (Figure [3]) we subtract flow from the 
measured phase speed to obtain the kink speed 
estimate (Figure 03). 

Now, since we have no information on the ex- 
pansion of the dark thread with height we make 
the initial assumption that the magnetic fleld is 
constant with height. This assumption is only for 
ease of calculation and provides an upper bound 
on the density scale height. As will be demon- 
strated in the next section, this assumption is not 
necessarily valid in such solar waveguides. 

First, we assume the density is gravitationally 
stratifled, such that 

pcxexp ^--^^^ , (8) 

where H^f / = H/ cos 6 is the effective density scale 
height due to the inclination of the d ark tube from 
the normal to the surface (see, e.g. lAschwanden 
I2OO4) . H is the density scale height and 6 is the 



inclination angle. We need the relationship 

where the subscript refers to the 'flrst' value of 
phase speed and density, Az is the difference in 
height between where Cko and Ck are measured. 
Substituting in the obtained values of phase speed 
for the second zero crossing with flow subtracted 
(Figure [2)3), where c^o is the phase speed at height 
km and is at height 7000 km, we obtain that 
Hef / = 2.4 ± 0.5 Mm corresponding to T 0.05 ± 
0.01 MK, where we have used the relationship 

H = ^^T, (10) 
prriHg 

and /cb is the Boltzmann constant, /i is the mean 
molecular mass, mn is the mass of hydrogen and 
g is the gravitational acceleration {2kB I J JimBQ ~ 
47 Mm/MK, see e.g., Aschwandeii] l2004l ). 

We can attempt to take into account expan- 
sion with height, which gives a new relationship 
between phase speed and scale- height, namely 

Using the values of expansion in the measured 
range for the larger UV/EUV jet, i.e. R/Ro ^ 
1.35 (from Figure [2^), and the same values of 
phase speed for the second zero crossing, we ob- 
tain Heff = 1.6 ± 0.5 Mm corresponding to T ~ 
0.03 ± 0.01 MK. Any larger value of expansion for 
the dark thread will give even smaller values of H 
and T, hence, assuming no magnetic fleld expan- 
sion along the thread provides an upper bound for 
density scale heights and temperature estimates. 
The obtained estimates for the temperature cor- 
respond to an intermediate temperature between 
chromospheric temperatures (10^^ K) and transi- 
tion region temperatures (10^ K). Co-spatial emis- 
sion with dark thread is observed in the 1600 A 
channel, which has contributions from transition 
region elements (i.e. C IV). Investigation of the 
emission reveals it appears to be transition region 
material associated with, rather than part of, the 
ejection of the dark thread. This suggests the ma- 
terial is at least cooler than the transition region 
plasma, i.e. T < 10^ K. 
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As we see signals of a reconnection event, this 
could suggest that the thread is a chromospheric 
loop that has been heated. The obtained estimates 
of temperature also agrees with the hypothesis 
given by our earlier arguments on why the filament 
appears dark in the 171 A channel, furthering cre- 
de nce to the s uggestion that we see a similar event 
to iLiu et al. I (jio09). Further, the obtained esti- 
mates for the density scale height are greater than 
the estimated scale heights derived for spicules 
from magnetoseismology in I Vert h et al. ( 2011 ) 



because the thread is hotter than Ca II H spicules. 

4.2. Estimation of flux tube expansion 

Now, the rest of the phase travel times occur 
during times of significant dynamic behavior in the 
larger jet (see. Figure [2]). We do not know if the 
dark thread also experiences similar dynamic be- 
havior. To make progress, we assume the thread is 
at least in a quasi-static state during the time for 
which the position of the maximum amplitude of 
the wave propagates through the cross-cuts. This 
means that the following magnetoseismological re- 
sults can be interpreted as a time average of the 
thread's structure. 

The phase travel time for the position of the 
maximum amplitude has a very similar fit when 
either a quadratic or exponential fit is used. We 
use an exponential fit since this is particulary con- 
venient for calculating the uncertainty in the rate 
of change of phase speed along the thread. This is 
necessary for quantifying the uncertainty in mag- 
netic field strength and plasma density gradients 
inferred from the following magnetoseismology. In 
Figure [7^ the measured phase speed is shown and 
in Figure [Tb we plot the measured phase speed 
minus the flow, i.e. what should be the actual 
phase speed of the wave, for the position of the 
maximum amplitude. 

Using Equation ([7]) it is obvious that we can 
write 

^0 V 

where and Ro are the values obtained at low- 
est height and and R are the values obtained 
at a greater height. Using the measured value of 
amplitude and the actual phase speed we can de- 
termine the the rate of flux tube expansion. Fig- 
ure [8^ shows the normalized expansion along the 
dark thread. 



What likely values of expansion should we ex- 
pect for the dark thread? It is possible, as men- 
tioned, that the dark thread will display the same 
behavior as the larger UV/EUV jet. The local ex- 
pansion with height of the dark thread is initially 
similar to that measured in the global jet struc- 
ture shown in Figure [2^. The large UV/EUV jet 
takes on the shape of a magnetic funnel, suggested 
as being the standard magnetic configuration for 
the quiet Sun (e.g. ICabriel I Il976l : IPeter I l2QQlh . 
Expected expansion factors from ext rapolations 
for qu iet Sun magnetic funnels, e.g. iTan et al. 
(|2Q1Q[ ). are 4-12 for heights less than 20 Mm 
and 1 — 3 for heights greater than 20 Mm. The 
calculated expansion in Figure [8^1 also appears to 
show the dark thread exhibiting a contraction of 
the cross-sectional area, with the profile being al- 
most symmetric about 3500 km. Is it possible 
this is a real contraction? Two likely scenarios 
suggest themselves if it is a real contraction: (i) 
if the thread is a chromospheric loop and recon- 
nection occurred close to a footpoint of the loop, 
then we may expect the dark thread to have a 
relatively symmetric profile about the apex; (ii) 
there are strong magnetic elements either side of 
the dark thread forcing the field lines to converge. 
However, the results from the previous section 
suggest the density is decreasing along the dark 
thread. It is possible this contraction has disap- 
peared (due to untwisting of the thread) between 
the time the maximum amplitude travels along the 
thread to the time the second minimum travels. 
The decreasing flow speed may be a sign that the 
thread has untwisted/opened up significantly be- 
tween these two events. However, the minimum 
amplitude also appears to show a decrease of am- 
plitude with height. 

Another option is that the kink wave experi- 
ences damping. Resonant absorption is known to 
be an efficient method for the dampings; of prop- 



agating kink waves (e.g., iTerradas et al. I I2OIOI : 
ICoossens et al. 2011) a nd is an ideal candidate for 
plasm a heating ( Erdelvill 19981 : iTarovan fc Erdelv~ 



|2009[ ). Also, dynamic behavior has been shown to 
attenuate/ amplify the amplitudes of MHD waves 
(iMorton Erdelyi I l2009bl: iMorton et al. I I2OIOI : 
Ruderman ll2010l : lErdelvi et al.ll201l[ ). This means 
the contraction would be an artifact as the os- 
cillation is damped. We can estimate the wave- 
length/ wave number of the kink wave. The typi- 
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Fig. 8. — (a) The estimated expansion of the flux tube (sohd hue) along the dark thread derived by means 
of magnetoseismology. The dashed lines correspond to the 95% confidence interval, (b) The estimated 
change of magnetic field strength (solid line) along the dark thread from magnetoseismology. The dashed 
fines correspond to the 95% confidence interval, (c) The estimated change in density (solid line) along the 
dark thread obtained by solar magnetoseismology. The short dashed line corresponds to the 95% confidence 
interval, whereas the long dashed line corresponds to a 70% confidence interval. 



cal period of the oscillations (e.g. see oscillation 
in Figure [5j) is 360 s. Using P = 27r/uj then, 
5400 < A < 18000 km and 5 x 10"^ < k < 
25 X 10-^ km-i for 15 < < 50 kms-^ This 
means 0.03 < kR < 0.15 for a radius of ~ 725 km 
(the radius is at resolution length scales) and res- 
onant damping is efficient at these wavelengths. 

We can also determine the magnetic and den- 
sity gradients along the dark thread. The mag- 
netic field varies as 

and this is shown in Figure [Hb. Further, using the 
values of the actual phase speed and the relation 

we can determine the density gradient. Figure [Hh. 
If we fit a decreasing exponential profile to the 
decreasing section of the density we obtain a scale 
height of H ~ 1.2 Mm and thread temperature of 
T ~ 2.6 X 10^ K. This temperature is less than 
the estimates obtained in the previous subsection 
and adds weight to the suggestion the thread is 
chromospheric in origin. Note, the reason for a 
70% confidence interval in Figure [8b rather than 
a 95% interval is due being unable to calculate 



higher confidence due to large relative errors on 
the density. 

In Figure [9] we provide a schematic representa- 
tion of the apparent scenario in the jet suggested 
by the magnetoseismology. It demonstrates the 
sub-resolution expansion of the thread and the 
density stratification that can be inferred from 
such an approach. 

5. Conclusions 

In the present study we demonstrate the suc- 
cessful application of solar magnetoseismology to 
a solar jet. In particular we determine informa- 
tion on various plasma parameters of a dark thread 
that is observed to occur with a UV/EUV solar jet. 
From the magnetoseismology we are able to esti- 
mate the temperature of the dark thread, the sub- 
resolution expansion and the magnetic and density 
gradients along the thread. It would be difficult 
to obtain this information through analysis of the 
SDO/AIA filter data or spectral lines profiles, e.g. 
Stokes profiles, line widths, etc. 

The nature of the dark thread is still a mys- 
tery as no similar feature has been reported pre- 
viously. The temperature estimates (T - 3 X 10^) 
obtained from the magnetoseismology suggest the 
plasma is cooler than the observable temperature 
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ranges of both the 304 A and 171 A. This is prob- 
ably the reason, along with the plasma being op- 
tically thick/dense, that the thread appears as a 
dark inclusion in the filters. A dense, T ~ 3 x 10^ 
plasma is likely to have originated in the chromo- 
sphere. The thread feature is possibly part of an 
emerging flux region that was responsible for the 
reconnect ion that caused the jet. The closest com- 
pa rison we find is the chro mospheric jet observed 
bv lLiu et al. 

the magnetoseismology, the 



(I2QQ9L 

With respect to 



fact we obtain reasonable temperature estimates 
for the plasma gives confidence that the density 
and magnetic field gradients also obtained are 
physical. The only issue that arises is the appar- 
ent contraction with height of the loop. A reason- 
able scenario that fits with the entire picture of 
the jet and the other results from the seismology 
is difficult to come by. It could be due to having 
neglected any influence of wave damping (e.g. due 
to resonant absorption) or because we are missing 
information on the dynamic evolution of the dark 
thread. Further theoretical studies into the influ- 
ence of dynamic plasma behavior on wave propa- 
gation are required if a full picture of an objects 
dynamic behavior is to be found via magnetoseis- 
mology alone. 
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Fig. 9. — Schematic representation of the sub- 
resolution features of the dark thread determined 
from magnetoseismology. 
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